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1 Introduction 

A magic square is an arrangement of a set of integers into a square array such that each row, column and the 
two diagonals add to the same value called the magic number. Magic squares generally fall into the realm 
of recreational mathematics [8|, |9| , however a few times in the past century and more recently, they have 
become the interest of more-serious mathematicians. In this paper we explore how permutation matrices 
play a key role in categorizing and enumerating magic squares. Also, we generalize the results of Mattingly 
[Tl] to find the determinant and eigenspectra of sets of magic squares that have various symmetries. 

Throughout the extensive history of magic squares (see [a]), the main interest lay in their construction. 
Only in the past 60 years has analysis of magic squares as matrices been performed. In 1917, Henry Ernest 
Dudeney [l| classified all of the 4 by 4 magic squares, mainly by symmetries. Figure [1] (shown on the 
following page) shows Dudeney's classification scheme. There are 12 classifications of magic squares of this 
size and the arcs shown in the figure for each type square signify the two cells that add to half the magic 
number. For example, the Diirer Magic Square, named after Albrecht Diirer (1471-1526), is 



(1) 



and it can be classified as a type III (or regular) magic square. This can be seen by adding up elements 
that are symmetric about the center or are linked in the Dudeney diagram. Each symmetric pair adds to 
17 (half the magic number of 34). 

In 1948, C. W. Trigg provided [I4] a way to calculate the determinant of each magic square in Dudeney's 
classification. In particular, he showed that each one of the 4 by 4 magic squares of Dudeney's type I- VI 
has determinant zero. In 2000, R. Bruce Mattingly, proved [7] that not only is the determinant of each 
type III or regular 4 by 4 magic square equal to zero, but this is true for any even-sized regular magic 
square. 

Mattingly's analysis used the reverse matrix to define a regular magic square, and used the definition 
to show some interesting properties about regular even-ordered magic squares. We expand the analysis to 
a larger set. In fact, we generalize Dudeney's type I- VI to higher-order magic squares and prove that they 
are all singular and show many properties of their eigenspectra. 
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2 Mathematical Background 

It is mathematically convenient to define magic squares in a more precise manner than discussed above. 
First, let the column vector e^ be the standard basis vector in M", the vector with all zeros except the ith 
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Figure 1: The classification scheme of Henry Ernest Dudeney [1|] for 4 by 4 natural magic squares. The 
arcs denote the pairs of cells that add up to half the magic number. In this case, the magic number is 
34. We will show that magic square of types I- VI, XI and XII satisfy equations containing permutation 
matrices and that generalized versions of types I- VI are singular. 



element, which is one. Also, let e be the vector in R"' with every element 1, and let J be the square matrix 
of size n defined by 



"l;J 



1, if i + j = n + 1, 
0, otherwise. 



(2) 



The matrix J, the identity matrix reflected either horizontally or vertically, is often called the reverse 
matrix since it is useful for reversing rows or columns of vectors and matrices. If A is a matrix, then JA 
is the matrix A reflected vertically and AJ is A reflected horizontally. A specific example of this is shown 
below. 

Since we are interested in understanding how permutation matrices transform magic squares, we define 
a magic square using matrix operations. 

Definition 1. A magic square of order n is an n by n matrix A consisting of real numbers with the 
following properties: 

(3) 
(4) 
(5) 
(6) 

where fi is called the magic number. 

If the numbers used in the magic square are the consecutive integers, 1,2,..., n^, then fi = n{n'^ + l)/2 
and the magic square is called natural. 



ef Ae = pt. 


fori = 1,2,.. 


■ ,n, 


e'^Aej = ^, 


for j = 1,2,. 


■ ,n, 


trace(A) = ^, 






trace(JA) = /i. 







Equations ([3]) and (JH) are the properties that the rows and columns respectively must sum to fi. 
Equations ([5]) and ([6]) are the properties that main and minor diagonals must also sum to fi. Writing the 
definition in this way allows proofs to proceed in a clearer manner. 

Although the numbers within magic squares have typically been defined on the integers (often it is 
important that they are non-repeating as will be indicated), the definition above and many of the results 
that follow are generally applied to magic squares on the reals. 

Since some of the results are applicable to semi-magic squares, we define that here. 

Definition 2. A semi-magic square is a n by n matrix A of reals that satisfies ^ and ^. If the 
elements of the matrix consist of a reordering of the numbers 1,2,... , n-^, then the result is called a natural 
semi-magic square. 

2.1 Symmetries of Matrices 

As is standard, a matrix A is symmetric if A = A-^. Permutation matrices, and in particular permutation 
matrices with particular symmetries, play an important role in this study. 

Definition 3. A square matrix A is persymmetric if it satifies J A = (JA) . A square matrix is 
bisymmetric if it is both symmetric and persymmetric. A square matrix is called singly symmetric if 
it is either symmetric or persymmetric but not both. 

The term persymmetric means that it is symmetric with respect to the minor diagonal. A result of 
these is that JAJ = A , which can be interpreted as the property that rotating the matrix halfway around 
is equal to its transpose. We will investigate bisymmetric permutation matrices in the next section. 

Another class of permutation matrices that play a role in magic squares are those with a 90° rotational 
symmetry. For example. 



(7) 



is such a matrix. More formally, we define the following: 

Definition 4. // a matrix A satisfies A'^ = J A, then A is said to be 90 "-symmetric. 

We will see later that permutation matrices that are bisymmetric or 90°-symmetric transform magic 
squares into other magic squares. 

2.2 Magic Squares formed by rotating and reflecting 

It is a standard known "fact" that for each magic square, there are an additional 7 magic squares found by 
reflecting and rotating the original one. Most who study magic squares consider the eight related squares 
as equivalent or isomorphic. However, if thought of as a matrix, then these 8 magic squares are different 
but still all magic squares. As has been seen, using matrix operations is important and the eight related 
magic squares can be written as A, JA, AJ, JAJ, A^, JA-^, A-^J, JA-^J. 
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For example, if A is the Diirer magic square in ([T]) and shown below, then 
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are the unique rotations and reflections of the Diirer magic square. Starting with any magic square of any 
order, it may seem obvious that these 7 transforms result in a magic square. We will show this below, not 
to prove an obvious result, but mainly to illustrate how the definition of the magic square in (l3])-(l6]) is 
useful in such proofs. 

Theorem 1. If A is a magic square then following are also magic squares: JA, AJ, JAJ, A , JA , 
A^J, JA^J. 



Proof. First, we show that 



5i) are satisfied for JA. Applying (^ results in 
ef (JA)e = (ef J)Ae = e^+i_,;Ae = //, 



for i = 1, 2, . . . , n. In this case efj = e^_|_^_^ because the reversed version of e^ is e„+i_ 
Similarly, 



e'^(JA)ej = (e'^J)Aej = e^^Aej 



fJ', 



for alH = 1, 2, . . . , n and e-^J = e-^. Property ([5]) applied to JA is ([6]) and vice versa. This shows that JA 
is a magic square and the proofs of the remaining matrices listed in the theorem are similar to this. 

D 



2.3 Eigenvalues and Eigenvectors of Magic Squares 

There are a few basics of eigenvalues and eigenvectors of magic squares. Much of this work appears in 
Mattingly [7] and is reproduced to demonstrate the proofs using the matrix definitions of the magic squares. 
Also, Loly et al. [6| gives an extensive review of the literature of both applicable matrix results as well as 
properties of magic squares. They emphasize eigenvalues and eigenvectors of general magic squares with 
a particular consideration on the characteristic polynomial. They also show many examples of specific 
low-order magic squares. 

Again, to illustrate the usefulness of the definition of the magic square given above, we show the 
following. 

Lemma 1. Let A be a magic square with magic number fi. 

• The vector e is an eigenvector of A with corresponding eigenvalue fi. 



The vector e^ is a left- eigenvector of A with corresponding eigenvalue /i. 



Proof. Left-multiply both sides of ([3]) by e to get eejAe = e/_i. This simplifies to Ae = fie, hence e is an 
eigenvector with associated eigenvalue /i. 

Similarly right-multiply both sides of <^ by e^ to get e-^Ae^e-^ = /xe^. This simplifies to e-^A = i^e^, 
hence e is a left-eigenvector with associated eigenvalue //. D 

Often, the magic number /i and vector e are called the magic eigenvalue and eigenvector. Mattingly [71] 
showed some results of eigenvalues and eigenvectors of regular magic squares and later we will generalize 
the results to other types of magic squares with certain symmetries. 

2.4 Even-Ordered Regular Magic Squares are Singular 

Trigg [ij] proved that all order-4, regular (or Dudeney type III) magic squares have zero determinanto 
Mattingly generalized this to any even-ordered magic square. We define a regular magic square and 
reproduce some of Mattingly's notation which is helpful later in this article. 

Definition 5. A regular magic square is an order-n magic square A with magic number fi that satisfies 

_ 2fx 



n 



,n. 



for every i = 1,2, ... ,n and j = 1,2, 

Mattingly [7| noted that a regular magic square A satisfies the property 



A + JAJ 
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-E, 
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where E = ee , the square matrix of all ones. This is an equivalent definition of a regular magic square, 
because the matrix JAJ is the matrix A refiected both vertically and horizontally and adding it to A 
results in a constant matrix. For example, if A is the Diirer magic square in ([T|, then 



A + JAJ 



In this case, /u = 34 and n = 4. 

Defining a regular magic square as a matrix A that satisfies dS]) allows a natural generalization of type 
III magic squares of Dudeney in Figure [1] to those of higher order and as we will show below in a similar 
manner to types I and II as well. 

Mattingly shows two main results. The first is that if A is an eigenvalue of a regular magic square, A 
with associated eigenvector u (and X ^ fj,), then —A is also an eigenvalue of A with associated eigenvector 
Ju. Secondly, he shows that a regular even-ordered magic square is singular. We will examine each of 
these theorems in more detail later in this paper as we generalize these results. 



^Although Trigg states that the sequence of integers used to create the magic square take on an arithmetic progression, his 
proof should extend to any non-natural magic square as defined here. 



3 The Role of Permutation Matrices 

As we have seen in both the definition of a magic square and a regular magic square, the reverse matrix J 
is an important permutation matrix in this context. Much of this paper explains transformations of magic 
squares based on permutation matrices, where J is just an example of one. 

In short, a permutation matrix that right (left) multiplies a matrix permutes the rows (columns) of a 
matrix. We saw above that if A is magic, then JA, AJ and JAJ is also magic. We are interested in what 
permutation matrices P result in PA, AP and PAP being magic. 

First, consider a square matrix A and a permutation matrix P of the same order. The matrix P can 
be considered as a map P where Pi is the column where the 1 appears in the zth row. The map can be 
considered a set of points of the locations of the Is. For example if 
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Then the map is P = {(1, 2), (2, 3), (3, 1), (4, 4)} or more compactly can be written (as in j3]) as (2 3 1 4). 
The matrix multiplication AP is a permutation of the columns of A under the map P. The matrix 
multiplication PA is a permutation of the rows of A under that map -P~^, which in this example is 
(3 1 2 4). Generally the matrix and the map can be used interchangeably, but typically we will denote 
when each is used. 

As we will see permutation matrices play an enormous role in magic squares, but first we will examine 
their role in semi-magic squares. 

Theorem 2. // P is a permutation matrix and A is a semi-magic square then PA, AP and PAP are 

semi-magic. 

Proof. We first show that PA is semi-magic by showing that ([3j) and ^ are satisfied. Let P be the map 
corresponding to the permutation matrix as defined above. Equation ([3]) applied to PA is 

ef (PA)e = (Pei)^Ae = e^Ae = ^, 

where k = P~ . This is satisfied for each k in {1, 2, . . . , n} since the map P is onto. Equation (j4]) applied 
to PA is 

e'^(PA)ej = (e'^P)Aej = e^Ae^ = fi. 

The proof that AP is semi-magic is similar, but not shown here, and the matrix PAP is semi-magic 
since both PA and AP are semi-magic. 

D 

Throughout this paper, we will use many examples of permutation matrices. In the order-4 case, 
there are 4! = 24 such matrices and when enumerated, they are done in a standard manner with Pi = 
(12 3 4), Ps = (1 2 4 3), Pg = (1 3 2 4), . . . , P22 = (4 2 3 1), P23 = (431 2), P24 = (432 1). Also, note 
that Pi = I, the identity matrix and P24 = J. 



3.1 Bisymmetric and 90°-Symmetric Permutation Matrices 

The permutation matrices that play the most significant role in transformation of magic squares are those 
that are bisymmetric. For example, I, J as well as 
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(9) 



are the only 6 order-4 bisymmetric permutation matrices. The subscripts are the standard enumeration of 
the permutation matrices. Here we study these matrices including both the count of bisymmetric matrices 
as well as showing their role in magic squares. 

Lemma 2. Let B{n) he the number of bisymmetric permutation matrices of order n. Then B[l) = 
1, B{2) = 2, S(4) = 6, B{n) = B{n - 1) if n > 3 is odd and B{n) = 2B{n - 2) + (n - 2)B{n - 4) i/ n > 6 
is even. 

Proof. The 1 by 1 matrix [1] is the only bisymmetric permutation matrix of this size and I and J are the 
only order-2 bisymmetric permutation matrices. In the order-4 case, I, J and the four matrices in Q are 
the only 6 bisymmetric permutation matrices. 

If n > 3 is odd, then in order to be bisymmetric, there must be a 1 in the center row, center column. Let 
P be such a matrix. Then the submatrix P([l, 2, . . . , (n — l)/2, (n + 3)/2, . . . , n)], [1, 2, . . . , (n — l)/2, (n + 
3)/2, . . . ,n)]) must also be bisymmetric, and this submatrix is order n — 1. Thus, the number of such 
matrices of odd order must be equal to the number of bisymmetric permutation matrices which is one 
smaller order or B{n — 1). 

Lastly, if n > 6 is even, then we construct all matrices of this order. If the 1 in the first row is 
in the first column (and thus by symmetry) also in the last row and last column, then the submatrix, 
P([2, . . . , n — 1], [2, . . . , n — 1]) must be bisymmetric and the number of these are B{n — 2). An additional 
B{n — 2) bisymmetric permutation matrices are found with a 1 in the first column last row (and again by 
symmetry) the last row and first column. 

The number of such matrices with the 1 on the first row and in column j with 2 < j < n — 1 is found 
by noting that this row is ej and by symmetry the first column is e^, the last column is e„+i_j and the 
last row if e^_,_^_ ■. The remaining submatrix to be filled is (n — 4) rows by (n — 4) columns and must 
be bisymmetric. Therefore the number of such matrices is (n — 2)B{n — 4) and adding to the number of 
matrices with a 1 in the upper left or upper right corner, the total number of even bisymmetric permutation 
matrices is B{n) = 2B{n - 2) + (n - 2)B{n - 4). D 

We will see later how the knowledge of the number of bisymmetric permutation matrices of a given 
order helps understand the count of the number of magic squares of a given order. We now prove a major 
result between general magic squares and permutation matrices. 

Theorem 3. If A is a magic square and P is a bisymmetric permutation matrix, then PAP is a magic 
square. 

Proof. Since P is bisymmetric, then P = P-^ and JP = (JP)^ or alternatively, JP = PJ, since J is 
symmetric as well. Also, from Theorem [21 PAP is semi-magic. Thus, we only need to prove the diagonal 
sum properties. Applying ^ and ^ to PAP results in 

trace(PAP) = trace(PPA) = trace(A) = fj, 
trace(JPAP) = trace(PJAP) = trace(PPJA) = trace(JA) = // 



7 



where JP = PJ is used, for any permutation matrix P ^ = P-'" and recalling that traces are invariant 
under cyclic permutations. (See [51]). 

D 

As an example let A be the Diirer Magic square in ([T|) and using the bisymmetric permutation matrix 
P3 given in (l9]), the result is 



P^AP, 



16 2 3 13 

9 7 6 12 

5 11 10 8 

4 14 15 1 



another magic square. Pre- and post-multiplying A by any of the other 5 bisymmetric permutation matrices 
results in other magic squares (or the same in the case of P = 1). 

Lemma 3. IfPis bisymmetric, then JP is bisymm,etric. 

Proof. We need to show that both (JP) = (JP)^ and (JJP) = (JJP)'^ or P = P^. These two statements 
comprise the definition of bisymmetric. D 

Above, we defined a 90°-symmetric matrix and gave an example of a 90°-symmetric permutation matrix 
in ([7]) and here we will show their role in transformations of magic squares and similar to that of bisymmetric 
permutation matrices, we will find the number of such matrices of a given order. 

Lemma 4. Let R{n) be the number of order-n permutation matrices which are 90 '-symmetric. Then 
R{1) = 1,R{2) = 0, -R(4) = 2, if n > 3 is odd then R{n) = R{n — 1) and if n is even and n > 6 then 
R{n) = {n-2)R{n-A). 

Proof. The only order-1 permutation matrix, [1] is also 90°-symmetric. There are two order-2 permutation 
matrices and neither are 90°-symmetric. The order-4 case can be shown by construction, let P be such a 
matrix. The first row cannot have a 1 in the first or last column, because this would imply that there is 
a 1 in each corner and such a permutation matrix does not exist. Thus, the first row can have a 1 in the 
2nd or 3rd column, and placing the 1 fills in the rest of the matrix. These two matrices are the matrix in 
([7]) and J times that matrix. 

Next, consider an odd order 90°-symmetric permutation matrix, P with n > 3. The 1 in the middle row 
of such a matrix must be in the middle column in order to be 90°-symmetric. The submatrix P[(l..(n — 
l)/2, (n -|- 3)/2..n), (l..(n — l)/2, (n -|- 3)/2..n), ] must also be a 90°-symmetric matrix, so R{n) = R{n — 1). 

Lastly, let P be an order-n 90°-symmetric permutation matrix with n even and n > 6. This proof 
continues by construction. As described above, the 1 in the first row cannot be in the first or last column, 
so a 1 can be placed in column 2 through n — 1 (for a total of n — 2 locations). Due to the the symmetry 
required, a one must also be located on the last column, last row and first column and zeros throughout 
the matrix in all rows and columns where the Is are located. This process fills in 4 rows and 4 columns 
symmetrically and the remaining submatrix (of size (n — 4) by (n — 4)) must be filled in with a 90°- 
symmetric permutation matrix. Thus the number of 90°-symmetric permutation matrices for an even 
order is R{n) = {n - 2)R{n - 4). D 

This lemma shows that the 90°-symmetric permutation matrices that exist are doubly-even order (size 
4/c) and odd order just larger (size Ak + 1). We also show that this proof leads to a construction of 90°- 
symmetric permutation matrices of a given size. As an example, let n = 8. We can select a 1 in first row 



(columns 2 through 7), we'h pick the 3rd column. Since the matrix need to be 90°-symmetric, this results 
in Is in locations to keep the matrix 90°-symmetric or 

10 0' 

00000001 


10 

.00000100. 

To fill in the remaining submatrix we can choose either the matrix in ([7]) or J times that matrix. Let's 
fill in the blank spots above with ([7]) to arrive at the 90°-symmetric permutation matrix: 

10 0' 
00010000 
1 
10 
10 
10 
10 
.00000100. 

The importance of these permutation matrices in the study of magic squares is due to the following 
theorem. 

Theorem 4. // P is a 90°- symmetric perm,utation matrix and A an order-n magic square, then PAP is 
a magic square. 

Proof. If P is 90°-symmetric, then it satisfies P^ = JP or alternatively, P = P-^J. From Theorem[2l PAP 
is semi-magic, so we need to only show the diagonal properties of the magic squares. Applying ([5]) to PAP 
we get, 

trace(PAP) = trace(PPA) = trace(PP^JA) = trace(JA) = // 

since traces are invariant under cyclic permutations and for any permutation matrix P^ = P^^. The last 
step arises from ([6]). In a similar manner, applying ([6]) to PAP is 



trace(JPAP) = trace(P'^AP) = trace(PP^A) = trace(A) 
where the last step is found in ([5]). 
Lemma 5. IfP is 90 "-symmetric, then JP is 90°- symmetric. 



^, 



D 



Proof. Applying the definition to JP leads to (JP) 
to JP = P^, the definition of 90°-symmetric. 



J(JP) or (JP)^ 



P. Taking the transpose leads 

D 



3.2 Extending the Family of Magic Squares 

The results of Theorems [3] and [J] indicate that we can extend the number of related magic squares beyond 
that given by Theorem [TJ From Theorem [1] and Lemmas [2] and U] there are 8{B(n) + R{n)) related magic 
matrices by applying PjAPj to the magic square A for Pj a bisymmetric or 90°-symmetric permutation 
matrix. However, as we will see double counting occurs and the number of unique matrices of the form 
PjAPj is only 4(S(n) + R{n)). 

To examine this more closely, let A be an order-4 magic square, and Pj is the ith standard permutation 
matrix. Theorem [3] indicates that lAI, P3AP3 PgAPg, P17AP17, P22AP22, and JAJ are also magic 
squares. In addition, PnAPn and P14AP14 are also magic from Theorem HI If the reflections and 
rotations in Theorem [1] are applied to these eight matrices, there are 64 matrices, howover only 32 are 
unique. For example, JAJ is listed in both Theorem [3] (since J is bisymmetric) and Theorem [TJ Another 
example shows that if one pre- and post-multiplies PnAPn by J (from Theorem[l]) the result is P14AP14, 
since JPn = P14. 

The double counting occurs because the permutation matrices listed above always come in pairs. If P 
is either bisymmetric or 90°-symmetric, then JP is as well, as was shown in Lemmas [3] and \E\ 

We will call the set of matrices starting with magic square A and formed by both reflections and rota- 
tions as well as transforming using the appropriate permutation matrix the family of transformations 
of A. For example, if we define A to be the Diirer Magic Square in ([T]) then A, PsAPsjPgAPg and 
PiiAPii or 



PsAP, 
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as well as the 8 rotations and reflections of each is the family of transformations of the Diirer magic square. 
In general, the following theorem explains this. 

Theorem 5. Let P be the set of bisymmetric or 90°-symmetric permutation matrices of order n, where only 
one ofP or JP is included in P. Let A be an order-n magic square arising from a set of non-repeating num- 
bers and Pj e P, then {PjAPj, JPjAPj, PjAPjJ, JPjAPiJ.PjA^Pj, JPiA^Pi,PiA^PiJ, JPjA^PjJ} 

is a set of unique magic squares. The size of this set is p{n) = 4(i?(n) + R{n)). 

Proof. For a given Pj, the 8 matrices listed is the matrix PjAPj and the 7 related matrices of Theorem [T] 
formed by reflecting and rotating. Each of these is unique since rotating and reflecting results in a unique 
matrix. Also, for a given A, the matrices PjAPj are unique, since each Pj is unique. Lastly, if Pj G P 
and JPj ^ P, then each rotation and reflection of PjAPj is unique. 

Since only half of the bisymmetric and 90°-symmetric permutation matrices Pj are included in the set 
P, the size of the set for a given order is half of the number of rotations and reflections (eight) times the 
number of permutation matrices in the set or B(n) + R{n), where these are given in Lemmas [2] and HI 

D 

The following table lists the number of unique magic squares in a family of transformations of a given 
magic square, for 3 < n < 10. Also, the total number of magic squares for a given order n will be a multiple 
of this number. 
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4(fi(n) + R{n)) 
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32 


32 


80 


80 


352 


352 


1248 



The total number of natural order-3 magic square matrices is 8, so the transformations account for all 
of them. In the case of order-4 magic squares, there are 7040 unique matrices, a factor of p(4) = 32, and 
Richard Schroeppel (given credit posthumously by Martin Gardner [2]) found that there are 2,202,441,792 
order-5 natural magic square matrices, also a factor of p(5) = 32. The number of higher-order magic 
squares are only estimates (see [lOt] for estimates for order 6), but p{n) gives a factor of magic squares 
and would provide a check on the accuracy and the family of magic squares in Theorem [5] can aid in the 
computation of all magics of a given order in that finding a single magic square results in finding all p{n) 
magic squares is its family. 

The factor of 32 for the order-5 magic squares was shown by Gardner {2] in that he states that in 
addition to the 8 reflections and rotations, there are two other transformations: 

1) "Exchange the left and right border columns, then exchange the top and bottom border rows." 

2) "Exchange rows 1 and 2 and rows 4 and 5. Then exchange columns 1 and 2, then columns 4 and 5." 

And thus he arrives at 4 • 2 • 2 • 2 = 32 isomorphic magic squares — four from the rotations, 2 from the 
reflections, and 2 each from the two above transformations. 

In light of our work here, the order-5 bisymmetric permutation matrices 



106 















1" 





1 

















1 

















1 





1 











0_ 



26 






1 








0] 


1 




















1 




















1 











1 


0_ 



where the numbers are the standard numbering system, are two of the 6 bisymmetric permutation matrices. 
If we transform an order-5 magic square A as P105 AP105 and P26AP26 the result will be the transforms 
1) and 2) listed above. According to Gardner, these two transformations can be composed, and the other 
factor of two (in the counting of isomorphic magic squares) results in applying Pio5P26AP26Pi05- The 
matrix multiplication 



P105P26 



is P45, a 90°-symmetric permutation matrix. Thus our counting of the order-5 isomorphic magic squares 
is equivalent to Gardner's. 
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3.3 Pandiagonal and Semipandiagonal Magic Squares 

There are two other categories of magic squares that are often studied. A pandiagonal magic square is 
one in which all of the broken diagonals also add to the magic number. For example. 



1 14 4 15 

8 11 5 10 

13 2 16 3 

12 7 9 6 



(10) 



11 



is a pandiagonal magic square since it is a magic square (with fi = 34) as well as aj^i + aj+2,2 + ^■j+3,3 + 
aj+4A = A* and a^j + aj+iA + Oj+2,3 + «j+3,2 = /^ for j = 1,2,3,4 and addition in the subscripts are 
performed modulo 4. (See [9| for more details about pandiagonal magic squares.) Examination of these 
matrices in the context of permutation matrices are found by the shifts of the identity matrix. That is, for 
order-4: 
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We can then use these permutation matrices to shift the broken diagonals to the main diagonals. For 
example, if A is the magic square in (llOp then 



AP 



10 



15 1 14 4 

10 8 11 5 

3 13 2 16 

6 12 7 9 



APioJ 



.4 14 1 15 

5 11 8 10 

16 2 13 3 

9 7 12 6 



The matrix APio has shifted the broken diagonal below the main diagonal to the main diagonal and APio J 
has shifted the broken diagonal above the minor diagonal to the main diagonal. When transformed as AP 
and APJ, the permutation matrices P17 and P19 shift the other broken diagonal to the main diagonal. 
In this light, we can formally define a pandiagonal magic square using permutation matrices. 

Definition 6. If A is a magic square with magic num,her fj,, then A is pandiagonal if A satisfies both 

trace(AP) = fi and trace(APJ) = //, 

for all P which are shifts of the identity matrix. 

Although the examples above used only order-4 matrices, the formal definition can be extended to any 
order. Pandiagonal magic squares of higher order are often studied. 

Another often-used classifying criteria is that of a semipandiagonal or semi-Nasik magic square. In this 
case only two of the broken diagonals are used — in the order 4 case only the one corresponding to P17 and 
P17J. Again, using permutation matrices, we can now extend the definition more formally. 

Definition 7. If A is a magic square with magic number /x, then A is semipandiagonal if 



trace(AP) 



^ 



and 



trace(APJ) = /i 



where P is the shift of the identity matrix with P symmetric. 



This definition allows the generalization of semi-pandiagonal magic squares to orders higher than 4. 
There is no permutation matrix of order 5 which is a shift of the identity matrix and symmetric. For order 
6, the permutation matrix that play the important role in the definition is: 

10 0' 
10 
1 
10 
10 
10 
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and an example of an order-6 semipandiagonal magic square is 
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where fi = 120. 

4 Classifying Dudeney Types using Permutation Matrices 

Next, we examine the work of Dudeney p| and his twelve order-4 types of magic squares which is shown 
in Figure [TJ We first investigate a simple question for all order-4 magic squares: if P is a permutation 
matrix, and Ag is a magic square of type q, for which permutation matrices is PA^P a magic square and 
what type is the result? In investigating this, we are able to generalize types I-VI to higher order magic 
squares and show properties of their determinants and eigenspectra. 

The directed graphs in Figure H] summarize the results of transforming all order-4 magic squares under 
PAP for a permutation matrix P. The upper-case letter types (A-D) listed above each graph are the same 
as the classification of Trigg [12 ] who grouped Dudeney's roman numeral types. Within each of the four 
graphs, the labelled arrows denote which sets of permutation matrices transform to other magic squares. 
The vertices of each graph are the Dudeney type. For example within type A, permutation matrices in the 
set A2 transform type I magic squares to type III magic squares (and vice versa). 

As can be seen in Figured! types I-III (type A) are always transformed within that set, as are types 
IV- VI (type B), VII-X (type C) and lastly the set within XI and XII (type D). As was shown in Theorems 
[3]and|4l the set of bisymmetric and 90°-symmetric permutation matrices transform magic squares to other 
magic squares. This set is either Ai (J ^44 or Ci (J C2 (J C3 |J C4. 

In the case of types I~VI, there are additional permutation matrices (the set A2 |J A3) that transform 
some magic squares to other squares. In addition, type VI magic squares do not transform in a consistent 
manner and thus to be fully understood has be subdivided into two sets VI', which are type VI magic 
squares that are semipandiagonal and VI", which do not have the property. 

4.1 Magic Classifying Permutation Matrices 

The transformations of type A and B magic squares can be generalized to higher-order magic squares. We 
will show that the generalized type A (and type B) magic squares can be transformed into other type A 
(and type B) magic squares using other permutation matrices. We first start defining the type I and II 
order-4 magic squares in a manner similar to that of Mattingly JTI] . If 
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(11) 



then if A/ is a type I magic square then A/ satisfies 



Ai + KA/K 



2/x 



n 



E, 



(12) 
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type A 



type B 



Ai,A2 




Ai,Ai 



Ai,Ai Ai,A^ 




Ai.Ai 



type C 



type D 





Ai (XI W 



A, 




^(XII) ) Ai 



Figure 2: A 
summary of the transformation of order-4 magic squares. If A is magic, and P is a permutation matrix, 
then the diagrams show for which permutation matrices is PAP magic and the Dudeney type of the original 
and transformed matrix. The vertices of the graphs are the Dudeney types and the arrows are labelled with 
the transformation sets: Ai = {Pi, Pg, P17, P24}, ^2 = {P2, P7, P18, P23}, ^3 = {Pe, Pio, Pi5, P19}, ^4 = 

{P3,Pll,Pl4,P22},Cl = {Pl,P24},C2 = {P3,P22},C3 = {P8,Pl7},C4 = {Pll,Pl4}, where the Sub- 

scripts are the index of the permutation matrices in standard form. The types (A-D) listed above each 
graph correspond to the types that Trigg [12] defined in terms of Dudeney's roman numerical types. Also, 
the set VI' is the type VI that are also semipandiagonal, whereas VI" is not semipandiagonal. 
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and if A// is a type II magic square then A// satisfies 

An + LA//L 



2^ 



E. 



n 



(13) 



Equations (J12p and ()13p are analogs of (|8]) for types I and II. We illustrate the relationship between the 
diagrams of Dudeney with this equation through an example. The matrix 

"1 14 4 15^ 

8 11 5 10 

13 2 16 3 

12 7 9 6 

is a natural order-4 magic square (/i = 34) that is type I. If we substitute this into (fT2|) . 



A + KAK 



1 14 4 15 

8 11 5 10 

13 2 16 3 

12 7 9 6 



+ 



16 3 13 2' 

9 6 12 7 

4 15 1 14 

5 10 8 11 



17 17 17 17" 

17 17 17 17 

17 17 17 17 

17 17 17 17 



^E, 
n 



where ^u = 34 and n = 4, so the type I matrix A satisfies (|12p . 

Magic squares of type IV, V, VI, XI and XII can also be written in terms of permutation matrices. Let 
Ag be a magic square of type q, then the following are true: 



A/y + LA/y = -^E, 
n 

Ay + KAy = ^E, 

n 
2/i 

Ay/ + JAy/ = — E, 

n 
2/i 
Axi + LAx/Ps = — E, 
n 

2/i 
Axil + KAx//P2 = — E. 
n 



(14) 
(15) 
(16) 
(17) 
(18) 



where P3 is defined in ([9|) and P2 = (1 2 4 3), written in compact notation. 

Initially, we examine types I-III which are written as equations using matrix operations in ([8]), (J12p 
and (J13p and note that K,L and J play an important role in these equations. Those matrices will be 
called magic classification permutation matrices in part because these permutation matrices can be used 
to classify a magic square in a way similar to that of Dudeney. 

Definition 8. A magic classifying permutation matrix or MCPM is a symmetric permutation 
matrix, P of order n and 

• if n is even then trace(P) = 0. 

• if n is odd then the middle element ofP is 1 and trace(P) = 1. 

The matrix J is the only MCPM of order 2 and 3 and the matrices J,L and K are the only magic 
classification permutation matrices of order 4. We note that K = P17 and L = Pg as shown in Q and 
the other two bisymmetric permutation matrices in Q, P3 and P22 are not MCPMs, nor are Pn (seen in 
([7]) and Pi4, the two 90°-symmetric permutation matrices. 

Before showing how to generalize the types I-III order-4 magic squares, it will be helpful to know how 
many MCPMs there are of a given order. 
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Lemma 6. Let C{n) denote the number of MCPMs of order n >2, then 



C{n) 



1)!! if n is even, 
2)!! if n is odd. 



Proof. We will prove this result by constructing even-order MCPMs, then showing the relationship between 
even and odd order. 

The only MCPM of order 2 is J, so C(2) = 1. Let P be a MCPM of order n > 4 and let j satisfy 
2 < J < n and be the column of the 1 in the first row of the matrix P; row 1 is e^ and column j is ei. 
(Note that to satisfy trace(P) = 0, the 1 cannot lie in column 1.) Also, by symmetry column 1 is ej and 
the jth row is e^. The remainder of the matrix is the submatrix P[(2..j — 1, j + l..n), (2..j — 1, j + l..n)] 
and must be a MCPM of order n — 2, which can be constructed in the same way. 

Thus the number of MCPMs of order n is the number of choices for placing the 1 in the first row times 
the number of MCPMs of order n — 2 or C{n) = (n — 1) • C{n — 2). Since C(2) = 1, for n an even number 
C{n) = {n-l)\\. 

If n is odd, the middle term must be a 1 by definition, and the submatrix resulting from removing the 
middle column and middle row is an order-(n — 1) MCPM. Thus the number of odd magic permutation 
matrices is the same as that of the even-ordered matrix just smaller than it or (n — 2)!!. D 

As an example, we show how to create a MCPM of order 6. We start with a 1 in the 4th column of 
the first row (and could have chosen any column between 2 and 6) and due to the fact that a MCPM is 
symmetric and a permutation matrix, so far this looks like: 

1 0' 








10 








The remaining submatrix consisting of the empty matrix elements above, must also be a MCPM. There 
are three choices for this: K, L, or J. To continue the example, consider L as the submatrix and the result 
is: 



1 0' 

10 

10 

10 

1 

10 



(19) 



This order-6 example is one of 5!! or 15 MCPMs. The first step in creating them is to place a 1 in one of 
the 2nd through 6th columns of the first row, thus 5 choices. By symmetry, there is a 1 in the same row 
of the first column. Then the resulting 4 by 4 submatrix can be either K, L or J. 



4.2 Classification of Magic Squares by Trigg Type 



Trigg [Ij] found the determinant of each of the order-4 880 magic squares (or 7040 unique matrices) based 
on categories of Dudeney types. In his proof, he classifies all 4 by 4 magics into four types (labelled A 
through D) as discussed above. Trigg's proof involves operations to change types II and III into a type I 
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magic square. Precisely, he says for a type II magic: Interchange 2nd and 3rd rows, and then interchange 
2nd and 3rd columns to get Type I. To show how this language is equivalent to using permutation matrices, 
recall that a type II magic square satisfies (J13p and let A// be a magic square of type II. Interchanging 
the 2nd and 3rd rows as well as the 2nd and 3rd columns of A//, a type II matrix, is A = P3A//P3 where 
P3 is defined in ([9]). When (fT3]) is multiplied both on the left and right by P3, the result is 

PsAj/Pg + P3LA,,LP3 = ^P3EP3 

n 

and using that P3L = KP3, LP3 = P3K as well as P3EP3 = E to get 

2/i 



P3 A77P3 + KP3 A77P3K 



-E. 



n 



This shows that the matrix A = P3A//P3 is a type I matrix from (J12p . 

The Trigg classification can be linked to the MCPM by the following definition. 

Definition 9. If A is a magic square that satisfies 

A + PAP = ^E 

n 



(20) 



for some MCPM P, then A is called a type A magic square. If the particular permutation matrix P is 
important, then it is called a type Ap magic square. 

In this language, the order-4 type I, II, and III magic squares can respectively be called type Ak, Al 
and Aj magic squares and any regular magic square is of type Aj . This definition also clearly allows the 
categorization of magic squares to that of any order. This and the next definition which generalizes type 
B magics is why we have named this set of permutation matrices, magic classifying permutation matrices. 

As an example of an order-6 type A magic square, we use a magic square (similar to that of Mattingly) 
consisting of the set {1,2,3,... ,38,39} in which the numbers 10, 20 and 30 have been removed. (A 
discussion of why we have done this is below.) The following magic square: 
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(21) 



is a type Ap, where P is the MCPM in ^ in that it satsifies A + PAP = (2^/n)E where fi = 120 and 
n = 6. 

Definition 10. If A is a magic square that satisfies either of the following: 



A + PA = ^E, 

n 

A + AP = ^E, 
n 



(22) 
(23) 



for some MCPM, P then A is called a type B magic square. If the specific permutation matrix is 
important, we call the matrix a type Bp magic square. 
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It is noted that types IV-VI satisfy equation (j22]) and on the surface (|23|) is unnecessary, however it 
is needed if magic squares are treated as matrices. Recall that Dudeney and Trigg assumed that a magic 
square that satisfies (j23p could be rotated to one that satisfies ()22p . In a similar manner to that mentioned 
above, the order-4 magic squares of types IV, V and VI can be called type Bl, Bk and Bj. 

Also, as an example of a type B order-6 magic square. 



(24) 



satisfies (I23p . where P is found in (I19p and /u = 120. This matrix is singular, has rank 4 and later we will 
show the rank and determinant of all type B magics. 

4.3 Type A Magic Squares 

In Theorems [3] and m we showed that if A is a magic square and P is either bisymmetric or 90°-symmetric 
that PAP is magic. This next theorem generalizes some of the results seen in the order-4 examples shown 
above. 

Theorem 6. // Q is a bisymmetric permutation matrix, A is a type Aj> magic square, and P' = QPQ 
then QAQ is a type A-pi magic square. 

Proof. The matrix QAQ is a magic square due to Theorem [3l Also, since A is a type Ap magic square, 
it satisfies (I20p . Multiplying (I20p on both the left and right by Q results in 

QAQ + QPAPQ = ^E 

n 

Letting P = QP'Q, and using QQ = I and QEQ = E, then 

QAQ + QQP'QAQPQQ = QAQ + P'QAQP' = ^E 

n 

which shows that QAQ is a type Ap/ magic square. 

Lastly, we need to show that P' is a MCPM, that is symmetric and has the same trace as P. First, 
P'^ = (QPQ)'^ = Q^P^Q^ = QPQ = P' and trace(P') = trace(QPQ) = trace(QQP) = trace(P) since 
traces are invariant under cyclic permutations. Thus P' is a MCPM. 

D 

In Figure HI we presented a set of directed graphs representing the transformations of order-4 magic 
squares by pre- and post-multiplication by permutation matrices. For example, a type I (or type Ak) 
matrix is transformed to a type II (or type Al) by P3. Using Theorem El we let P = K in (llip and 
Q = P3 in ([9]), then P' = QPQ = P3KP3 = L, therefore using P3, the type Ak (type I) matrix is 
transformed to a type Al (type II) matrix. 

There are permutation matrices other than bisymmetric ones that transform type A magic squares to 
other type A magic squares. For example, P2 (which is symmetric, but not bisymmetric) transforms a 
type Ak (type I) to a Aj (type III) — see Figure HI From the results in Figure HI the set of permutation 
matrices that transform type A magic squares have the property that the permutation matrix is singly 
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symmetric. However, moving to higher-order magic squares, this is not true. It appears from examining 
type A magic squares with sizes above 4 by 4, that further analysis is needed. 

The following result shows that any type Ap magic square can be transformed to any other type Ap/ 
magic square. 

Theorem 7. /f P and P' are even-ordered MCPMs, then there exists a symmetric permutation matrix Q, 
such that P' = QPQ. 

Proof. Let the size of the matrices be 2n and V and V' be the graphs that correspond to the permutation 
matrices P and P'. Since P and P' are even-ordered MCPMs, the graphs V and V' consist of n edges, each 
which connects a vertex of degree 1. These two graphs are isomorphic and from Goodaire and Parmenter 
[J], there exists a permutation matrix Q such that P = QP'Q-^, where Q is the permutation matrix that 
relabels the vertices of V to that of V' . Because of the structure of these graphs, the relabeling is done in 
a symmetric way — that is, the relabeling occurs as swaps. Thus the permutation matrix Q is symmetric 
and the theorem holds. 

D 

This theorem gives further insight into the categorization and properties of magic squares. For example, 
since there exists a matrix that converts a type Ap magic square to a type Ap/ magic square, this implies 
that there are equal numbers of magic squares of Pj-type for every MCPM, Pj. This is because for each 
type Ap magic square, there is a transformation to a type Ap/ magic square. In the order-4 case, there 
are 384 matrices of each type I, II and III (of Dudeney's classification). 

Theorem [7] also shows how to select Q to convert from a type Ap magic square to a type Ap/ magic 
square. Consider an example where 
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(25) 



We build the following graphs V (on the left) and V' (on the right): 





The relabeling of the first graph as swapping vertices 2 and 4 results in the second graph. Thus the 



19 



permutation matrix 



Q 



1 0' 
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10 
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10 
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(26) 



perforins the relabeling and this matrix satisfies P = QP'Q for P and P' in (j25p . There are also 3 
additional swaps and permutation matrices that will also convert P to P', thus Q is not unique. 
The matrix 
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(with fi = 120) is type Ap, where P is given in 



A' = QAQ 



Using Theorem El then the matrix 

32 16 4 27 39 2^ 

5 25 37 19 22 12 

23 7 31 14 11 34 

13 38 1 8 36 24 

26 6 29 17 9 33 

21 28 18 35 3 15 

where Q is the permutation matrix in (j26p . is a type Ap/ magic square for P' defined in (j25p . 

5 Singly-Even Ordered Magic Squares 

Many of the examples in this article use either order-4 magic squares, mainly due to the fact that there 
are only 7040 unique matrices, and such a small number allows analysis of all magic squares of this order. 
Less work has been done for order-5 magic squares, but it is known that odd-ordered magic squares have 
significantly different properties from those of even-ordered ones, and as we will discuss in this section 
singly-even (those of order 4/c + 2 for k G Z+) magic squares have different properties from those of 
doubly-even (order 4A; for k € Z^) magic squares. 

Also, it is noted that most of the order-6 magic squares that we have shown as examples here are 
not natural ones — instead we have selected squares that arise from removing 10, 20 and 30 from the set 
{1, 2, 3, ... , 38, 39}. Examples include those in (j2ip and (|24p . The main reason for selecting a nonnatural 
magic square is the main result in this section, that no natural type A magic square exists. 

The main result below relies on a proof of Planck [111], who showed that a natural singly-even pandi- 
agonal magic squares does not exist. A consequence of this, as he pointed out, is that there are no regular 
magic squares. As noted in the next theorem, this can be extended to any type A natural magic square. 

Theorem 8. There are no type A natural magic squares of singly-even order. 
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Proof. Assume that A is a type Ap natural magic square of singly-even order. From Theorem [7] there is a 
matrix Q such that P = QJQ. Then from Theorem[6l QAQ is a type Aj (or regular) magic square, which 
is a contradiction, since Planck [ill] proved that no such natural regular magic square exists. Therefore no 
type A natural magic square exists. 

D 

6 Determinants, Eigenvalues and Eigenvectors of Even- Ordered Magic 
Squares 

We now explore determinants, eigenvalues and eigenvectors of type A and B magic squares. In particular 
we will show that magic squares of these types are singular. 



6.1 Type A Magic Squares 

In this section, we investigate type A magic squares, defined above. The theorems presented here are 
similar to those in Mattingly [7|] but are generalized from a regular magic square to any type A magic. 
First, we define. 



A-^E. 
n 



(27) 



The matrix Z is used in the proofs of theorems below, however is also useful for exploring the symmetry 
of a magic square. For example, let 

20 42 30 41 33 52 3 39' 



23 45 26 13 62 24 32 35 
16 7 19 50 56 40 60 12 
47 28 48 10 59 38 



22 8 
63 64 



4 14 29 31 44 11 



57 43 27 55 6 37 17 18 
1 2 54 34 21 51 36 61 

58 49 53 25 5 15 9 46 



(28) 



which is a type A magic square. It satisfies A + PAP = (2/i/n)E where 



10 0' 
10 
00000001 
10 
10 
00010000 
10 
10 



(29) 



and fj, = 260. This permutation matrix can be constructed in a manner similar to that discussed in Lemma 
Eland this is one of 7!! or 105 order-8 MCPMs. 
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Figure 3: A Dudeney diagram for a magic square of type Ap, where P is found in (|29p . The magic square 
A in (I28p has a symmetry in which the entries connected by the arcs in the diagram add to 2;u/n = 65 in 
this example. 



Using (l27|l and ()28[) . this example results in 
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15.5 


-22.5 


26.5 


5.5 


30.5 


31.5 


-28.5 


-18.5 


-3.5 


-1.5 


11.5 


-21.5 


24.5 


10.5 


-5.5 


22.5 


-26.5 


4.5 


-15.5 


-14.5 


-31.5 


-30.5 


21.5 


1.5 


-11.5 


18.5 


3.5 


28.5 


25.5 


16.5 


20.5 


-7.5 


-27.5 


-17.5 


-23.5 


13.5 



The matrix Z can be used to develop a Dudeney-like diagram for the matrix P in (j29p . which is fascinating 
in that there is no matrix P in the definition of Z. To create such a diagram, one draws links between 
numbers of the same magnitude and opposite signs — the diagram is found in Figure [3l 
This next theorem shows the eigenvalues of Z are related to those of A. 



Theorem 9. Let A be a type A even-ordered magic square and let Z he the matrix defined in (21). Then 
A and Z have the same eigenvalues, except that /x is replaced by in the spectrum ofZ. Also, both Z and 
A have the same eigenvectors. 

The proof of this theorem can be found in Mattingly [7|. His proof applied only to regular magic 
squares, however with minor changes, this extends to any type A magic square. 

Theorem 10. Let A be a type Ap magic square. If X ^ n is a eigenvalue of A associated with eigenvector 
X, then —X is also an eigenvalue with associated eigenvector Px. 

— Z. Left- and right- 
(30) 



Proof. Let Z be the matrix defined in (I27p . First, we will show that PZP 
multiplication of (p7|) by P results in, 



PZP = PAP 



n 



PEP, 
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and using ((20]) as well as PEP = E. Adding (|27]) and ((30]) and using PEP = E results in 

Z + PZP = A - -E + PAP - -E. (31) 

n n 



Using ([20]), the right hand side of 1^ is 0, so PZP = -Z. Since P = P"\ Z is similar to -Z, thus Z 
and — Z have the same eigenvalues. A consequence of this is that if A 7^ is an eigenvalue of Z, then —A 
is also an eigenvalue of Z. 

Next, assume that x is an eigenvector of Z with eigenvalue A. Multiply PZP = — Z on the right by x 
to get 

PZPx = -Zx = -Ax 

and multiplying this through by P^^ = P, the result is 

ZPx = P-I(-Ax) = -APx, 

so —A is an eigenvalue of Z with related eigenvector Px. Lastly since Z and A have the same eigenvalues 
(except for fi) and eigenvectors as seen in Theorem [9l then the statement of the theorem is true. D 

In light of Theorems [9] and [TU\ we examine A in (I28p in more detail. The eigenvaluecl are 

{260, 0, -61.80, 61.80, 40.17i, -40.17i, 11.39 i, -11.39 i} 

and interesting properties include: 

• The number 260 is the magic eigenvalue as is stated in Lemma [1] 

• The matrix is singular since is an eigenvalue. We will prove this about all type A magic squares 
below. 

• Each of the other eigenvalues arise as pairs of equal magnitude and opposite sign. This is a conse- 
quence of Theorem [TOl 

If we extend the example shown above to include eigenvectors, consider the eigenvectoio associated 
with A4 = 61.80 which is 

X4 = (0.230, 0.717, -0.355, -0.342, 0.237, -0.020, -0.344, -0.123)"^ 

and the eigenvector associated with A3 = — A4 = —61.80 is 

X3 = (0.717, 0.230, -0.123, -0.020, -0.344, -0.342, 0.237, -0.355)"^ 



and it can be seen that X3 = PX4 with P defined in (|29p . recalling that left-multiplying by P interchanges 
pairs of elements in a vector. The remaining pairs of eigenvalues and eigenvectors are related in a similar 
manner. 

Theorem 11. Suppose that a square matrix B is similar to — B. // a Jordan Block J^ (A) appears m times 
in the Jordan Canonical Form ofB then so does J^ (—A). 



^The eigevalues in this example and later related eigenvectors are calculated numerically. The results presented are rounded 
to either 3 or 4 significant digits to illustrate the theorems. 

^The eigenvectors shown here are found numerically as was the eigenvalues and are selected such that the 2-norm equals 1 
with the same signs. 
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Again, this theorem comes from Mattingly [71] and is true because both B and — B have the same 
Jordan Form. The following now extends Mattingly's proof that regular magic squares are singular to any 
type A magic square. 

Theorem 12. If A is a type A magic square of even-order, then det(A) = 0. 

Proof. Let Z be defined in (|27|) . For each nonzero eigenvalue A of Z, Theorem 1111 ensures that the Jordan 
blocks of A are paired with those of —A. 

Then the sum of the algebraic multiplicities of the nonzero eigenvalues of Z must even. Note that 
is eigenvalue of Z, so its multiplicity must also be even. If Z has at least two eigenvalues that are 0, this 
implies at least one eigenvalue of A is 0. Thus, det(A) = 0. D 



6.2 Type B Magic Squares 

Theorem 13. // A is an order-2n magic square that satisfies either 
Furthermore, the rank of A is at most n + 1 . 



or /123\) then A is singular. 



Proof. We first consider a matrix A that satisfies (I22p and let Pi be the column where the one appears 
in the ith row of the matrix P in (j22p . There are n pairs {i,Pi) with i < P,. Thus the n row operations 
Ri + Rp^ -^ Rp^ results in n rows contain the row vector (2/i/n)e^. 

Since there are n identical rows, the determinant is 0. Also row operations can be used to get n — 1 
rows of zeros, thus the dimension of the null space is at least n — 1 and thus the rank of the matrix can be 
no larger than n + 1. 

If A satisfies ()23p then the same argument can be made using column operations to generate n columns 
containing the vector (/i/n)e. Similarly the rank of these matrices is also at most n + 1. 

D 

We demonstrate this result with a natural 8 by 8 magic square. The following matrix 



14 39 63 16 33 24 40 31' 



56 


47 


42 


19 


15 60 


1 


20 


44 


54 


10 


7 


53 22 


8 


62 


59 


17 


4 


35 


52 37 


29 


27 



6 48 61 30 13 28 36 38 ' ^ ' 

21 11 55 58 12 43 57 3 

9 18 23 46 50 5 64 45 

.51 26 2 49 32 41 25 34 

is a type B magic — that satisfies (|22p — with P = J. The magic number in this case is // = 260. Applying 

the row operations Ri + R% ^ R%^ R2 + Rj ^ Rj .^ R'i + Rq ^ Re and i24 + -R5 — > R5 leads to 

14 39 63 16 33 24 40 31' 

56 47 42 19 15 60 1 20 

44 54 10 7 53 22 8 62 

59 17 4 35 52 37 29 27 

65 65 65 65 65 65 65 65 

65 65 65 65 65 65 65 65 

65 65 65 65 65 65 65 65 

65 65 65 65 65 65 65 65 
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If we also perform the row operations —R5 + R^ ^ Re, —R5 + Rj ^ Rj and —R^ + R% ^f Rs, the last 3 
rows are zeroed out. This shows that the rank of the matrix is at most 5. 

Also, the eigenspectrum of 1^ is {260,0,0,0,-53.8553,49.6710,2.0921 ± 43.6941i}. Recall that the 
number of zero eigenvalues of a matrix is equal to the dimension of the null space of the matrix and for an 
order-n square matrix A, dim(A) + rank(A) = n, so the rank of this matrix is 5. 

7 Magic Square classification and eigenspectra 

Above we discussed that the rotations and reflections of a magic square as well as left- and right-multiplying 
by a bisymmetric and 90°-symmetric permutation matrix generates a family of magic squares, however only 
the bisymmetric ones apply to the next theorem. 

Theorem 14. Let {Pj} be the set of hisymmetric permutation matrices that is the same size as a magic 
square A. Then 

• A, A-^,PjAPj, PjA-^Pj each have the same eigenspectrum for each i. 

• JPiAPj,PjAPjJ,PjA PjJ, JPjA Pj each has the same eigenspectrum for each i. 

Proof. Recall that A and A have the identical eigenspectrum {A}. 

Since Pj is a symmetric permutation matrix, it satifies Pj = P~ so PjAPj is similar to A and thus 
has the same eigenvalues. PjA^Pj is also similar to A-^, so also has the same eigenvalues. 

We now show the latter four matrices listed in the theorem have the same eigenspectrum. Let the 
eigenspectrum of JPjAPj be {A}. Since J satisfies J~^ = J, then left- and right-multiplying JPjPPj by 
J is PjPPjJ, therefore this has the eigenspectrum {A}. 

Lastly, the transpose of JPjAPj is 

(JP.APi)^ = Pf A^Pf J^ = P.A^PiJ 

and it has the same eigenspectrum of JPjAPj. A similar analysis shows JPjA^Pj has the same eigen- 
spectrum. 

D 

This theorem proves that half of the magic squares in a set of this form of magic squares have one set 
of eigenvalues and the other half has a (possibly) different set. 

8 Conclusions 

The seminal work of Dudeney on categorization of order-4 magic squares and the subsequent work of Trigg 
on the determinant of each order-4 magic square can be put into a greater framework. We have extended 
the family of 8 magic squares formed by rotation and reflection to include transformations by classes of 
permutation matrices. In addition, many of the types of order-4 magic squares that Dudeney described 
can be naturally extended to higher order magic squares. 

As has been noted in many works on magic squares, there are fundamental differences between even 
and odd-ordered magic squares and there are further differences between singly- and doubly-even magic 
squares. We extended the proof that there are no natural singly-even regular magic squares to the fact 
that there are no natural singly-even type A magic squares. 

In addition, we generalized the work of Mattingly, who show that all even-ordered regular magic squares 
are singular. We extended this to all even-ordered types A and B magic squares as well showing some 
properties of eigenvalues and eigenvectors of type A magic squares. 
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